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ABSTRACT 

The broad-line radio galaxy 3C 111 has been suggested as the counterpart of the 7-ray source 
3EG J0416+3650. While 3C111 meets most of the criteria for a high-probability identification, 
like a bright flat-spectrum radio core and a blazar-like broadband SED, in the Third EGRET 
Catalog, the large positional offset of about 1.5° put 3C111 outside the 99% probability re- 
gion for 3EG J0416+3650, making this association questionable. We present a re-analysis of all 
available archival data for 3C111 from the EGRET archives, resulting in detection of variable 
hard- spectrum high-energy 7-ray emission above 1000 Me V from a position close to the nominal 
position of 3C 111, in three separate viewing periods (VPs), at a 3a level in each. A second vari- 
able hard-spectrum source is present nearby. At >100 Me V, one variable soft-spectrum source 
seems to account for most of the EGRET-detected emission of 3EG J0416+3650. A follow-up 
Swift UVOT/XRT observation reveals one moderately bright X-ray source in the error box of 
3EG J0416+3650, but because of the large EGRET position uncertainty, it is not certain that the 
X-ray and 7-ray sources are associated. Another Swift observation near the second (unidentified) 
hard 7-ray source detected no X-ray source nearby. 

Subject headings: gamma rays: observations — galaxies: active — galaxies: individual(3C 111) 


1. Introduction 

One of the main scientific goals of the upcoming 
7-ray astronomy satellite mission GLAST (Ritz 
2007; Michelson 2007) is to shed light on the 
nature of powerful relativistic extragalactic jets, 
which are ejected from the nuclei of some active 
galaxies (AGN). Originally, this class of AGN was 
defined based on the bright and prominent ra- 
dio emission from these jets: the radio-loud pop- 
ulation of active galaxies. Based on data from 
EGRET, the high-energy 7-ray telescope on the 
Compton Gamma Ray Observatory , it was realized 
(Fichtel 1994; Thompson et al. 1994) that the 
largest population of extragalactic ~GeV 7-ray 
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sources is represented by those radio-loud AGN 
whose jets are pointed at a small angle to the 
line of sight. The fact that the bright compact 
radio emission of these so-called blazars makes 
excellent targets for parsec-scale resolution very- 
long-baseline interferometric (VLBI) observations 
of their jet structure ties an intimate link between 
radio VLBI and high-energy 7-ray astronomy. 

Of special interest are the questions of where 
in the AGN jets the bright 7-ray emission is pro- 
duced, and how the emission is interacting with 
its immediate environment and with other parts 
of the jet. This knowledge would enable us to put 
crucial constrains on the processes of jet forma- 
tion, collimation, and acceleration. GLAST is ex- 
pected to yield densely sampled 7-ray light curves 
of hundreds of extragalactic jets that are bright 
enough to be detected on time scales of days to 
weeks, and thousands on time scales of months 
to years. Most of these objects will be blazars; 
in fact, all but two of the firmly identified extra- 
galactic EGRET sources (LMC, Sreekumar et al. 
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1992; CentaurusA, Sreekumar et al. 1999) are 
blazars. In contrast to those, radio galaxies have 
larger inclination angles and the beaming is there- 
fore weaker. That allows better (deprojected) lin- 
ear resolution with VLBI observations and in the 
case of stratified jet structures, it allows observa- 
tions of the slower jet layers, e.g., a sheath, whose 
emission may be swamped by the much brighter 
beamed emission from faster jet regions, e.g. a 
fast spine, in blazars. Ghisellini et al. (2005) have 
presented such a spine-sheath stratified-jet model, 
predicting detectable 7-ray emission from the nu- 
clei of some radio galaxies. In their model, each of 
the two emission regions sees photons coming from 
the other part relativistically enhanced because of 
the relative speed difference, giving rise to an ad- 
ditional inverse-Compton emission component. 

In this letter, we report on the identification 
of the broad-line radio galaxy 3C 1 1 1 as a 7-ray 
source, responsible for a portion of the source 
3EG J0416+3650. Such an association was sug- 
gested as possible in the third EGRET catalog 
(Hartman et al. 1999, 3EG), but was considered 
unlikely because of the large positional offset of 
3C111 from 3EG 0416+3650. The present work 
was stimulated by the recent report of Sguera 
et al. (2005), who used multi wavelength data to 
strengthen the case for the association between 
3C 111 and 3EG 0416+3650. 

Using detailed analysis of reprocessed EGRET 
data (Esposito et al. 1999; Bertsch et al. 2001), 
we show here that 3EG J0416+3650 is most likely 
composed of at least two, and more likely three, 
separate sources, one of which is in good positional 
agreement with 3C 1 1 1 . As demonstrated recently, 
based on VLBA monitoring data at A2cm from 
the MOJAVE program (Kadler et al. 2008), the 
parsec-scale jet of 3C 111 shows a variety of phys- 
ically different regions in a relativistic extragalac- 
tic jet, such as a compact core, super luminal jet 
components, recollimation shocks and regions of 
interaction between the jet and its surrounding 
medium, which are all possible sites of 7-ray pro- 
duction. Its relatively large inclination angle of 
~ 19° makes 3C 111 a particularly well-suited tar- 
get for tests of structured-jet models, such as the 
model of Ghisellini et al. (2005). We describe our 
re-analysis of the available EGRET data on 3C 111 
in Sect. 2 and discuss our results and their impli- 
cations in Sect. 3. 


2. Analysis 

2.1. EGRET Data Analysis 

We have investigated the possibility that 
3EG J0416+3650 is a superposition of at least 
two sources, one being at or near the position 
of 3C111. The analysis for 3EG was based on 
photons with energies >100 MeV, for which the 
68% containment angle of the spectrum-dependent 
point spread function (PSF) is always several de- 
grees. For energies >1000 MeV, the PSF is con- 
siderably smaller, slightly less than one degree. 
On the other hand, the number of photons with 
energies above 1000 MeV is much smaller than 
that above 100 MeV. Both of those integral en- 
ergy ranges were investigated for this study. 

For each viewing period (VP) in which 3C111 
and 3EG J0416+3650 were in the field of view 
(FoV) of EGRET, the following procedure was ini- 
tially used: 

1. Using the standard EGRET likelihood soft- 
ware (like), two sources were modeled, one 
at the 3EG position for 3EG J0416+3650 
and one at the nominal position of 3C 1 1 1 , 
and like then provided the best estimates of 
the significance and flux of the two assumed 
sources. 

2. If either of the two assumed sources had a 
test statistic (TS) of at least 2.5 (roughly 
1.6a), like was then allowed to optimize 
the significance of the two sources by moving 
their positions. 

The results from this procedure are shown in 
Table 1, for the two energy ranges and the two 
assumed sources. The first five columns show the 
VPs, the starting dates of the observations, the 
lengths of the VPs, EGRET wide (W) or narrow 
(n) FoV mode, and the off-axis angle of 3C111. 
The line showing VP=P1234 is the analysis of 
the summed data for the VPs included in 3EG, 
VP0002 through VP4270. The EGRET team has 
normally used data out to 19° in narrow-field 
mode, but as indicated in Column 5, that limit was 
stretched for the present study. This means that 
we were using regions of very low exposure, with 
possible slight systematic distortion of positions 
(Esposito et al. 1999). For each source and en- 
ergy range, the table shows the TS at the starting 
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position, the change in position (from the starting 
point) to obtain the maximum TS, and the result- 
ing photon flux in units of 10“ 8 photons cm -2 s -1 . 
Blank entries for TS indicate values of 0.00. For 
TS values less than 2.5 (roughly 1.6cr significance), 
there are no entries in the next two columns. 

The density of detections in this table is 
very low, even with the low requirement of 
TS>2.5. Three of the 3C111 detections (P1234 
and VP4270 for E>100MeV and VP0150 for 
E>1000MeV) have offset position values larger 
than would be expected, and out of line with 
the other detections. The VP0002 detection 
(>100 MeV) of J0416 is probably suspect for the 
same reason. 

Only one VP (3211) provides a reasonably 
strong detection of 3EG J0416+3650, which agrees 
well with the information provided in 3EG. The 
summed exposure PI 234 provides about the same 
significance of detection, indicating that the ob- 
ject was emitting 7 rays at a low level even when it 
was not significantly detectable in individual VPs. 
The two weak detections of 3EG J0416+3650 at 
E>1000MeV (VP4270 and P1234) agree well in 
position with the E> 100 MeV, but show very low 
photon fluxes. 

After excluding the suspect detections, there 
are six VPs that support the proposed 3C111 
7-ray source, three in E> 100 MeV and three in 
E> 1000 MeV. In each of these six cases, the value 
of offset in position is compatible (considering the 
energy range) with the position of 3C111. Only 
one VP (0310) shows detections in both energy 
ranges; however, each of the detections in that 
VP is the strongest in its energy band. Like 
3EG J0416+3650, 3C 111 appears to radiate 7 rays 
with a low duty cycle, occasionally becoming fairly 
bright. Note that the best exposure, VP0390, does 
not show detection of either of the two assumed 
sources. 

The number of photons >100 MeV from 3C 111 
was 64, 44 of them in VP0310. Only about 14 pho- 
tons >1000 MeV were detected from 3C 111, six in 
VP0150, four in VP0310, and four in VP3250. For 
3EG J0416+3650, there was a total of 152 pho- 
tons >100 MeV, 40 of them in VP3211, the most 
significant single detection. PI 234 yielded seven 
photons above 1000 MeV from 3EG J0416+3650, 
but as indicated in Table 1, they produced a TS 
of only 3.43, about 1.8cr. 


We have further investigated a small region 
around 3EG J0416+3650 by generating 5° x 5° 
likelihood maps, for the energy ranges >100 MeV 
and >1000 MeV, for each VP in which 3C 111 and 
3EG J0416+3650 were in the FoV. A sample of 
the more informative likelihood maps is shown in 
Figure 1. Maps for the summed data, P1234, are 
shown in Figure 2. We draw the following conclu- 
sions from these maps: 

1. There appear to be at least three sources 
(possibly more) contributing, two in the 
>1000 MeV energy range and at least one 
in the >100 MeV range. One of the two 
>1000 MeV sources is very near the position 
of 3C 111; for the other, which is of com- 
parable strength, we find no obvious iden- 
tification (see below). In the following, we 
designate this source EGEV J0426+3747. 

2. The two >1000 MeV sources appear to have 
very hard spectra, since there is little or no 
significant emission detected in the lower en- 
ergy bands. 

3. The >100 MeV emission is dominated by a 
source near the 3EG catalog position. This 
source seems to have a very soft spectrum, 
since there is no significant detection in the 
>1000 MeV energy range. In the following, 
the name 3EG J0416+3650 refers implicitly 
to this soft source. 

4. The duty cycle of all three sources seems to 
be in the range 10-25%. 

5. The >1000 MeV map for VP0150 shows why 
the position offset for 3C111 is so large 
(0.423°) in Table 1 in that VP. There are 
only about seven photons near 3C111 in 
that map. One of them is nearly a degree 
away from 3C 111. This shifts the maximum 
likelihood position toward that single pho- 
ton, which could be either from the tail of 
the PSF, or from the diffuse Galactic fore- 
ground. Thus this VP, which had previously 
been ignored, provides additional support 
for the >1000 MeV detection of 3C 111. 

As an additional check on the reality of the two 
>1000 MeV sources, we have compared the num- 
ber of occurrences of >2 a excesses in 17 maps 
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Table 1 

EGRET observations of 3C 111 and 3EG 0416+3650 analyzed in this study. 


>100 MeV > 1000 MeV 

3C 111 J 0416+3650 3C 111 J0416+3650 


VP 

Start Date 

Days 

FoV 

OfT-Axis a 

[deg] 

TS C 

Offset 

[deg] 

Flux b 

TS d 

Offset 

[deg] 

Flux b 

TS C 

Offset 

[deg] 

Flux b 

TS d 

Offset 

[deg] 

0002 

04-22-1991 

6 

W 

24.82 


_ 


12.14 

1.196 

50 


_ 


0.88 

_ 

0005 

05-04-1991 

3 

W 

22.82 

0.15 

- 


1.21 

- 



- 



- 

0010 

05-16-1991 

14 

W 

29.31 


- 



- 



- 



- 

0150 

11-28-1991 

14 

W 

10.01 

0.04 

- 


2.91 

0.299 

8 

8.96 

0.423 

2.0 

0.19 

- 

0310 

06-11-1992 

14 

W 

20.79 

10.01 

. 0.19 

26 

0.55 

- 


8.86 

0.188 

3.2 


- 

0365 

08-12-1992 

8 

W 

6.44 

1.64 

- 


0.73 

- 



- 



- 

0390 

09-01-1992 

16 

W 

5.45 


- 


1.58 

- 


1.30 

- 


0.03 

- 

2130 

03-23-1993 

6 

W 

20.72 


- 


0.42 

- 



- 



- 

2210 

05-13-1993 

11 

W 

25.79 

0.66 

- 



- 



- 



- 

3211 

02-08-1994 

7 

W 

20.60 


- 


16.90 

0.185 

55 


- 


2.41 

- 

3215 

02-15-1994 

2 

W 

20.60 

6.76 

0.45 

64 


- 



- 



- 

3250 

04-26-1994 

14 

W 

14.47 

1.71 

- 



- 


4.14 

0.092 

2.3 

0.78 

- 

4120 

02-28-1995 

7 

W 

25.41 

0.16 

- 



- 



- 


0.25 

- 

4260 

08-08-1995 

14 

w 

22.82 


- 


0.94 

- 



- 



- 

4270 

08-22-1995 

16 

w 

7.91 

2.98 

1.85 

20 

0.13 

- 


0.20 

- 


3.81 

0.049 

P1234 

- 

- 

w 

— 

2.73 

1.22 

4 

16.48 

0.064 

10 

8.59 

0.109 

0.8 

3.43 

0.031 

5020 

10-17-1995 

14 

n 

28.65 


- 


0.74 

- 



- 



- 

7289 

10-13-1998 

21 

n 

23.80 


- 


- 

- 



- 



- 

8290 

09-14-1999 

14 

n 

21.13 

4.06 

0.08 

77 


- 



- 



- 

9185 

04-25-2000 

14 

W 

25.65 

0.09 

- 



- 


0.31 

- 



- 


a Off-axis angle for the position of 3C 111 
b Flux in units of 10 -8 ph cm -2 s _1 
c Test statistic at the position of 3C 111 
d Test statistic at the position of 3EG J0416+3650 


4 



(5° by 5°) with the statistically expected value. 
The results are shown in Table 2. For the re- 
gions of the maps away from the two ” source” po- 
sitions, the number of >2 a excesses is about as 
expected from statistical fluctuations, taking into 
account the number of PSFs contained within each 
map. Near the locations of the two >1000MeV 
sources, however, there are far more excesses than 
expected, which provides added confidence in the 
detections of 3C111 and EGEV J0426+3747. 

2.2. Swift Follow-Up X-Ray Observations 
and Their Analysis 

There are no known high-energy sources close 
to the nominal best-fit positions of 3EG J0416+3650 
and EGEV J0426+3747. Therefore, we have con- 
ducted follow-up optical/UV and X-ray observa- 
tions with the UVOT and XRT instruments on 
the NASA satellite Swift (Gehrels et al. 2004). 
The 3EG J0416+3650 field was observed with 
XRT for 10.4 ksec on December 14, 2007 2 , and 
the EGEV J0426+3747 field for 9.0ksec on 20 De- 
cember, 2007 3 . No significant X-ray source was 
detected in the 23.6 x 23.6 arcsec XRT FoV around 
EGEV 0426+3747 above the background count 
rate level of 0.00045 ctss -1 . The XRT observa- 
tion of the 3EG J0416+3650 field detected one 
4.8a source at the position 04 h 15 m 54.3 s, +36 d 
49m 26 s, with an uncertainty of 5.8 arcsec. De- 
spite the fact that this position is very close to the 
best position estimate for 3EG 0416+3650, the 
large position uncertainty of the EGRET source 
precludes a firm association between them. 

We extracted spectra from the source (which 
we designate Swift J041554.3+364926) and back- 
ground regions using XSELECT, and response and 
ancillary files provided by the Swift calibration 
data base (CALDB) 4 . The low count rate of 
0.0026 ± 0.0006cts/s -1 limits the ability to de- 
rive spectral information. However, it is striking 
that only 4 photons (out of 26) have been de- 
tected between 0.2 keV and IkeV. To characterize 
the spectrum, we used Cash statistics (Cash 1979) 
in XSPEC and determined the maximum allowable 
power-law photon index in a power-law fit under 
the assumption of zero absorption. The result is 


2 Observation ID 00031057001 

3 Observation ID 00031072001 

4 http : //heasarc . gsf c . nasa.gov/docs/caldb/swift/ 


slightly dependent on the data binning chosen, 
but suggests a very hard spectrum, with the pho- 
ton index in the range 1.1 -1.3. If absorption is 
allowed, the photon index is poorly constrained, 
and steeper values up to 2.7 are statistically al- 
lowed. 

2.3. Spectral energy distributions 

The non-contemporaneous spectral energy dis- 
tribution (SED) of 3C111 is shown in Figure 3. 
The millimeter flux range was determined (with 
permission) from the SMA Submillimeter Calibra- 
tor List archived data 5 (Gurwell et al. 2007). The 
R-band flux range has been determined from data 
obtained in a long-term monitoring program with 
the Liverpool telescope (PI: I. McHardy) between 
November 2004 and July 2007 (S. Jorstad, priv. 
communication) . 

The 3C 111 SED is similar to those of EGRET- 
detected flat- spectrum radio quasars (FSRQs), ex- 
cept that the rollover in the NIR, optical, and UV 
is steeper here. This makes the ’’valley” between 
the synchrotron and Compton ’’bumps” appear 
deeper. This could be due to absorption, but it 
must be intrinsic to 3C111, since the data points 
have been corrected for local absorption. 

Figure 4 shows (a) a UVOT image of the 
Swift J041554. 3+364926 region, with a circle con- 
taining the XRT X-ray photons superimposed; 
and (b) an SED containing the very limited data 
available. The EGRET points shown are for com- 
parison only; despite its proximity to the 3EG 
position for the 7-ray source, it is not certain 
that Swift J041554. 3+364926 is actually the same 
object as 3EG J0416+3650, because the EGRET 
position uncertainty is quite large (about 0.63° 
in radius). Furthermore, 3EG J0416+3650 is 
quite variable, and in the SDSS POSS-I Blue 
Survey image from a 1955 Palomar plate, there 
is an additional object present about 25 arcsec 
west of Swift J041554. 3+364926 (at 04:15:52.411; 
+36:49:24.77). This object does not appear in any 
other optical or IR image we have located. 

In the Swift J041554.3+364926SED, the UVOT 
data are plotted twice, with and without dered- 
dening, because it is not certain whether the ob- 
ject is nearby or extragalactic. The open circles 


5 


5 http : //smal . sma . hawaii . edu/callist/callist . html 



Table 2 

EXCESSES >2(7 IN 17 MAPS (5° BY 5°) OF THE FIELD AROUND 3EG J0416+3650. 


VP 

Full 

3C 111 

0426+3747 

Background 


Field a 

Field b 

Field" 

Field d 

0002 

0 

0 

0 

0 

0005 

s 

0 

s 

0 

0010 

1 

0 

0 

1 

0150 

2 

1 

1 

0 

0310 

1 

1 

0 

0 

0360 

1 

0 

0 

1 

0365 

0 

0 

0 

0 

0390 

1 

0 

0 

1 

2130 

1 

0 

0 

1 

2210 

1 

0 

1 

1 

3211 

2 

0 

1 

1 

3215 

s 

0 

s 

0 

3250 

4 

1 

1 

2 

4120 

2 

0 

0 

2 

4260 

0 

0 

0 

0 

4270 

1 

0 

0 

1 

8290 

1 

0 

1 

0 

Sum 

18 + 2s 

3 

5 + 2s 

11 

Rate 

10% 

18% 

29% 

7% 


Note. — s (in VPs 0005, 3215, columns 2 and 4) - a 
single photon >1000 MeV very near the source position, 
significant because of the short exposures; a 187 PSFs; b 
17 PSFs; c 17 PSFs; d 153 PSFs 


6 


Table 3 

3C111 HISTORICAL SED DATA 


Frequency* 1 

Band 

Instrument 

Flux b 

Flux Error 

Reference 

(0.24 - 24.1) x 10 23 

> 100 MeV 

EGRET(max) 

1.92 x 10“ 10 

7.07 X 10 -11 

this work 

(2.41 - 24.1) x 10 23 

> 1000 MeV 

EGRET(max) 

1.11 x 10“ 10 

7.00 X 10 -11 

” 

(2.41 - 24.1) x 10 23 

> 1000 MeV 

EGRET (avge) 

2.90 x 10 -11 

1.50 X 10 -11 


(0.24 - 24.1) x 10 23 

> 100 MeV 

EGRET(avge) 

3.00 x 10 -11 

1.90 X 10 -11 


(0.24- 24.1) x 10 23 

> 100 MeV 

EGRET(VP3215) 

4.28 x 10“ 10 

2.06 X 10“ 10 


(2.41 - 24.1) X 10 23 

> 1000 MeV 

EGRET(VP3215) 

< 1.53 x 10“ 10 


r> 

(2.41 - 7.23) X 10 20 

1-3 MeV 

COMPTEL 

< 1.30 x 10“ 10 


Maisack et al. (1997) 

(7.23 - 24.1) x 10 20 

3-10 MeV 

COMPTEL 

< 1.30 x 10“ 10 


” 

(2.41 - 7.23) X 10 21 

10 - 30 MeV 

COMPTEL 

< 3.00 x 10“ 10 



(1.81 - 2.41) X 10 20 

0.75- 1 MeV 

COMPTEL 

< 2.90 x 10“ 10 



(1.20 -3.60) x 10 19 

50- 150 keV 

OSSE 

3.52 x 10 -11 

6.56 x 10“ 12 

Johnson et al. (1997) 

(9.70 - 24.0) x 10 18 

40 - 100 keV 

ISGRI 

8.10 x 10 -11 

7.40 x 10“ 12 

Beckmann et al. (2006) 

(4.80 - 9.70) X 10 18 

20 - 40keV 

ISGRI 

6.27 x 10 -11 

5.70 x 10“ 12 


(3.39 - 47.2) x 10 18 

14- 195 keV 

BAT 

1.25 x 10“ 10 

9.00 x 10“ 12 

Tueller et al. (2007) 

(4.80 - 19.3) x 10 18 

2 — 80 ke V 

BeppoSax 

4.99 x 10" 11 

9.98 x 10” 12 

Kadler (2005); Kadler et al. (2008) 

(4.80 - 24.0) X 10 17 

2-10 keV 

BeppoSax 

2.49 x 10 -11 

4.98 x 10” 12 


(1.20 - 4.80) X 10 17 

0.5 — 2 keV 

BeppoSax 

4.36 x 10“ 12 

8.72 x 10“ 13 

” 

(4.80 - 24.0) x 10 17 

2- lOkeV 

XMM 

6.57 x 10 -11 

3.94 x 10“ 13 

” 

(1.20 - 4.80) x 10 17 

0.5 — 2 keV 

XMM 

1.38 x 10 -11 

3.94 x 10“ 13 

” 

(4.80 - 24.0) x 10 17 

2 - lOkeV 

ASCA 

3.92 x 10 -11 

7.84 x 10“ 12 


(1.20 - 4.80) x 10 17 

0.5 — 2 keV 

ASCA 

7.42 x 10“ 12 

1.48 x 10“ 12 

” 

1.14 x 10 15 

UVW1 

UVOT 

< 6.98 x 10“ 13 


F. Verrecchia, priv. comm. 

8.57 x 10 14 

U 

UVOT 

< 6.95 x 10“ 13 



6.92 x 10 14 

B 

UVOT 

1.00 x 10“ 12 

(+1.0, -0.5) x 10“ 12 

” 

5.55 x 10 14 

V 

UVOT 

2.30 x 10“ 12 

(+1.1, -0.7) x 10“ 13 

” 

4.29 x 10 14 

R 

Liverpool 

(4.22 - 7.25) x 10“ 12 


S. Jorstad, priv. comm. 

2.40 x 10 14 

J 

2MASS 

1.19 x 10 -11 

4.00 x 10“ 13 

2MASS archive Skrutskie et al. (2006) 

1.82 x 10 14 

H 

2MASS 

1.59 x 10 -11 

4.62 x 10“ 13 

” 

1.36 x 10 14 

K 

2MASS 

1.92 x 10 -11 

3.66 x 10“ 13 

» 

2.50 x 10 13 

12 micron 

IRAS 

2.43 x 10 -11 

3.65 x 10“ 12 

Golombek et al. (1988) 

1.20 x 10 13 

25 micron 

IRAS 

2.69 x 10 -11 

4.04 x 10“ 12 


5.00 x 10 12 

60 micron 

IRAS 

1.61 x 10 -11 

2.42 x 10“ 12 

” 

3.40 x 10 11 

0.85 mm 

SMA 

(28.4- 5.27) x 10“ 12 


SMA archive 1 

2.30 x 10 11 

1.3 mm 

SMA 

(25.9 - 0.94) x 10“ 12 



9.00 x 10 10 

3 mm 

SMA 

(3.66 - 0.91) x 10“ 12 



1.50 x 10 10 

2 cm 

UMRAO 

(1.05-0.30) x 10“ 12 


Kadler et al. (2008) 

8.00 x 10° 9 

3.6cm 

UMRAO 

(6.40 - 3.20) x 10“ 13 



5.00 x 10° 9 

6 cm 

UMRAO 

(4.50 - 3.00) x 10“ 13 


)5 


a Frequency in units of Hz 
b Flux in units of erg cm -2 s -1 
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Table 4 

3EG J 0416+3650 SED data 


Frequency 51 

Band 

Instrument 

Flux b 

Flux Error 

Reference 

(0.24 - 24.1) x 10 23 

> 100 MeV 

EGRET (avge) 

7.03 x 10 -11 

1.89 x 

io - 11 

this work 

(0.24 - 24.1) X 10 24 

> 100 MeV 

EGRET(VP3211) 

4.42 x 10 -10 

1.29 x 

IQ- 10 

” 

(3.39 - 47.2) x 10 18 

14- 195 keV 

BAT 

< 1.78 x 10 -11 



Tueller et al. (2007) 

(4.82 - 24.1) X 10 17 

2 - lOkeV 

XRT 

1.05 x 10- 13 

5.24 x 

IO" 14 

this work 

1.48 x 10 15 

UVW2 

UVOT(observed) 

4.43 x 10~ 14 

1.48 x 

io - 14 

55 

1.34 x 10 15 

UVM2 

U V OT (observed) 

5.38 x 10“ 14 

1.35 x 

10~ 14 

55 

1.14 x 10 15 

UVW1 

UVOT (observed) 

1.94 x IO -13 

2.28 x 

10 -14 

” 

8.57 x 10 14 

U 

UVOT(observed) 

7.80 x 10" 13 

3.43 x 

10" 14 

55 

6.93 x 10 14 

B 

U V OT (observed) 

9.29 x IO -13 

5.57 x 

10 -14 

55 

5.55 x 10 14 

V 

UVOT(observed) 

1.13 x 10- 12 

9.40 x 

IQ- 14 

” 

1.48 x 10 15 

UVW2 

UVOT(dereddened) 

5.16 x 10~ 12 

1.72 x 

10 -12 

55 

1.34 x 10 15 

UVM2 

UVOT(dereddened) 

8.06 x 10“ 12 

2.02 x 

IO -12 

55 

1.14 x 10 15 

UVW1 

U VOT(dereddened) 

6.25 x 10~ 12 

7.38 x 

IQ- 13 

55 

8.57 x 10 14 

U 

UVOT(dereddened) 

1.06 x 10 -11 

4.67 x 

10 -13 

55 

6.93 x 10 14 

B 

UVOT(dereddened) 

8.15 x 10- 12 

4.89 x 

10 -13 

55 

5.55 x 10 14 

V 

UVOT(dereddened) 

6.22 x 10 -12 

5.17 x 

10" 13 

55 

4.29 x 10 14 

R 

2MASS 

< 4.58 x 10~ 12 



2MASS archive Skrutskie et al. (2006) 

2.40 x 10 14 

J 

2MASS 

2.66 x 10“ 12 

2.66 x 

io - 13 

” 

1.82 x 10 14 

H 

2MASS 

1.73 x 10“ 12 

2.60 x 

io - 13 

55 

1.36 x 10 14 

K 

2MASS 

1.57 x 10 -12 

1.57 x 

10- 13 

55 

1.40 x 10 09 

20 cm 

NVSS 

< 1.40 x 10~ 17 



Condon et al. (1998) 


a Frequency in units of Hz 
b Flux in units of erg cm - 2 s -1 


denote the observed fluxes; the dereddened fluxes 
are shown as filled circles. Note that the fall-off 
in the optical-UV is extremely steep, especially 
in the observed fluxes. If Swift J041554. 3+364926 
is an extragalactic source, this can be attributed 
to absorption in a fairly dense foreground molec- 
ular cloud, part of the Taurus complex, with ab- 
sorbing column density Nr — 2.18 x 10 21 cm -2 
(Kalberla et al. 2005). NED shows Ay of 1.893m 
for the catalog position of the EGRET source, 
about 0.1° away. In the extragalactic scenario, 
Swift J041554.3+364926 must be a radio quiet 
source because of the very low upper limit on its 
radio flux density, 0.7 mJy at 1.4 GHz, from the 
NVSS survey (Condon et al. 1998). The shape of 
the dereddened UV/optical/IR range suggests an 
AGN ’’blue bump”, usually attributed to thermal 
emission from an accretion disk. However, with- 
out far IR and/or submillimeter data points, it is 
not possible to be certain of this interpretation. 
Based on previous experience, a radio-quiet AGN 
is an unlikely identification for an EGRET source. 

Alternatively, if Swift J041554. 3+364926 is 
within our Galaxy, it would need to be an M star 


or a cataclysmic variable (CV) in order to pro- 
duce the observed X-rays (R. Osten and S. Drake, 
personal communication). In order to be visible 
optically, an M star would necessarily be in front 
of the molecular cloud, which is 140 pc distant 
(Elias 1978), so the observed IJVOT data points 
should be considered, rather than the dereddened 
ones. The observed X-ray emission is stronger 
than expected from an M star except during an 
exceptional flare (R. Osten and S. Drake, personal 
communication). If the object is a CV, its lumi- 
nosity is much greater, so it must be behind the 
Taurus absorber, making the dereddened UVOT 
points appropriate; in that case the X-ray emis- 
sion is reasonable. Based on previous experience, 
neither an M star or a CV would be a likely iden- 
tification for an EGRET source. 

Table 3 and Table 4 provide the values plotted 
in Figures 3 and 4, with references. 

3. Conclusions 

The 3EG source 3EG J0416+3650 seems to be 
composed of at least three highly variable sources, 
two hard-spectrum (3C 111 and EGEV J0426+3747) 
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and one soft-spectrum. One of the hard-spectrum 
sources is close to the radio galaxy 3C 1 1 1 , 
and presumably associated with it; the other, 
EGEV J0426+3747, has no obvious identification 
at other wavelengths. The soft-spectrum gamma- 
ray source may be associated with the XRT source 
Swift J041554.3+364926. The region under study 
is 8 — 10° from the Galactic plane, so one or both 
of these additional sources could be local. 

Our re-analysis of the available EGRET data 
supports the previous association of 3EG J041 6+3650 
with the broad line radio galaxy 3C111, but 
with 3C111 responsible for only a portion of 
the 7-ray emission. Furthermore, it explains the 
relatively large positional offset noted in Hart- 
man et al. (1999). We have compiled a his- 
torical SED of 3C111 which shows that the X- 
ray data may well extrapolate into the EGRET 
range, particularly during flares, which is in agree- 
ment with the intermittent nature of detections 
in the individual EGRET viewing periods. It 
is interesting to note that we detect 3C111 at 
almost exactly the gamma-ray flux that is pre- 
dicted by equation (12) in Ghisellini et al. (2005): 

(1.41 14.1) x 10 -11 ergs -1 cm -2 , scaling from 

the 5 GHz values in Table 3. Note that 3C111 
(z=0.0485) is the most distant radio galaxy de- 
tected in 7 rays, about twice as far as NGC 6251 
(Mukherjee, Halpern, & Mirabal; z=0.0247). 

The detection of strong and variable gamma- 
ray emission from 3C 111 further supports the hy- 
pothesis that radio galaxies may represent an im- 
portant class of GLAST sources. GLAST’s LAT 
detector is about a factor of 30 more sensitive 
than EGRET, with a factor of ~3 better PSF 
above 1 GeV. It is scheduled for launch in mid- 
May, 2008, and will continuously scan the en- 
tire sky, so it will be far more sensitive to tran- 
sient events than EGRET. In the case of 3C 111, 
EGEV J0426+3747, and Swift J041554.3+364926, 
the smaller PSF of the LAT will clearly separate 
these three sources and allow us to measure their 
time- variable gamma-ray spectra and duty cycles. 

We thank Goro Sato, Francesco Verrecchia, 
Svetlana Jorstad, and Chris Shrader for their 
help and advice in compiling the SED data for 
this study. Steve Drake and Rachel Osten have 
kindly provided insight regarding X-ray emis- 
sion from stars. MK has been supported by 


the NASA Postdoctoral Program at the Goddard 
Space Flight Center, administered by Oak Ridge 
Associated Universities through a contract with 
NASA. This research has made use of data ob- 
tained from the High Energy Astrophysics Science 
Archive Research Center (HE AS ARC), provided 
by NASA’s Goddard Space Flight Center, and 
the NASA/IPAC Extragalactic Database (NED), 
operated by the Jet Propulsion Laboratory, Cal- 
ifornia Institute of Technology, under contract 
with NASA. We have also made use of data prod- 
ucts from the Two Micron All Sky Survey, a joint 
project of the University of Massachusetts and the 
Infrared Processing and Analysis Center, Califor- 
nia Institute of Technology, funded by NASA and 
NSF. Submillimeter flux ranges for 3C111 were 
obtained from the on-line archives of the Submil- 
limeter Array (SMA), a joint project between the 
Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and As- 
trophysics, funded by the Smithsonian Institution 
and the Academia Sinica. 
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Fig. 1. — Likelihood maps for the seven viewing 
periods with TS>3 (anywhere in the map, either 
energy range). The left column panels are for 
E>100MeV, while the right column panels show 
E>1000MeV, with fewer photon but considerably 
smaller point spread function. Crosses show the 
positions of 3C 111 and 3EG J 0416+3650; ref. Fig- 
ure 2. 


10 



Galactic Latitude 



164 163 162 161 160 

Galactic Longitude 


Fig. 2. — Likelihood maps for the summed expo- 
sure of all viewing periods, for >100 MeV (top) 
and >1000 MeV (bottom). 



iog(frequency, Hz) 


Fig. 3. — Spectral energy distribution for 3C111. 
The data points are non-contemporaneous; their 
values and references are tabulated in Table 3. 
Thick vertical lines indicate historical ranges for 
radio and submillimeter observations. The thick 
vertical line in the optical indicates the range of R- 
band fluxes from 2004-2007 observations with the 
Liverpool Telescope - see text. The NIR through 
UV points have been dereddened; see text. For 
the EGRET data, the average value are indicated 
with open diamonds, the most significant single 
viewing period, VP0310, uses downward-pointing 
open triangles, and the maximum flux, VP3215, is 
shown with upward- pointing triangles. The lowest 
upper limit for a single VP (not shown) is similar 
to the average. 
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Fig. 4. — (a) UVOT image (UW1 filter) of the 
region around Swift J041554. 3+364926. The cir- 
cle indicates the position uncertainty of the X- 
ray source; (b) Spectral energy distribution for 
Swift J041554. 3+364926. The data points are non- 
contemporaneous; their values and references are 
tabulated in Table 3. For the UVOT data points, 
the observed values are shown as open circles, 
while dereddened values are shown as filled cir- 
cles. The EGRET values may not be associated 
with the lower- frequency source - see text. 
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